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Abstract 

We compute the wall velocity in the MSSM. We therefore generalize the SM equa- 
tions of motion for bubble walls moving through a hot plasma at the electroweak 
phase transition and calculate the friction terms which describe the viscosity of 
the plasma. We give the general expressions and apply them to a simple model 
where stops, tops and W bosons contribute to the friction. In a wide range of 
parameters including those which fulfil the requirements of baryogenesis we find 
a wall velocity of order v w ~ (5 — 10) • 10~ 2 much below the SM value. 
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1 Introduction 



The generation of the baryon asymmetry of the Universe in a first order phase 
transition (PT) of the electroweak theory is attractive because all of its ingredients 
could be tested in high energy experiments in the near future. However, it is for 
sure now that there is no strong first order PT in the Standard Model (SM) 
with a Higgs mass above the experimental lower bound, indeed there is no PT 
at all [1]. The SM is very successful, but there is common agreement that it 
has to be extended to a more general theory. That this theory will contain 
supersymmetry is still controversial but in lack of alternatives this is a very useful 
hypothesis. In the MSSM and also in an extension with an additional singlet field 
(NMSSM) [2-4] it is possible to obtain a strong first order PT without violating 
the experimental Higgs mass bounds. In the MSSM this is related to a scalar 
partner of the right handed top which is very light in the symmetric phase [5-7] . It 
increases the cubic term in the effective 1-loop scalar potential. In the NMSSM 
a SHiH2 term arises already at tree level and acts effectively as a 3 -term if 
the vacuum expectation value (vev) <S> is comparable to the Higgs vevs [8]. 
These types of models have also more freedom to realize CP-violating effects. It 
is well known that supersymmetric models allow for spontaneous CP-violation 
at T = 0. The parameter space is strongly restricted by experimental bounds 
on the electric dipole moment of the neutron (for a discussion and references 
see e.g. [9]). Moreover to produce a baryon asymmetry at the electroweak scale 
CP-violation within the bubble wall is needed. Therefore it has been proposed 
that a temperature induced transitional CP- violation might occur [10-12]. In 
[13, 14] it was shown that in the MSSM spontaneous CP-violation does not occur 
throughout. Even with maximal explicit CP violating phases the variation of the 
corresponding phase in the Higgs system is strongly suppressed. 

The baryon asymmetry arises in a two step process: first the expanding wall 
sweeps through the hot plasma separating Higgs phase and symmetric phase 
with a CP-violating spatially varying Higgs vev. It generates an asymmetry 
between left handed quarks and their antiparticles diffusing in front of the bubble 
starting from an asymmetry between stops, tops, charginos, neutralinos and their 
antiparticles. In a second step this asymmetry in the hot plasma in front of the 
bubble wall is transformed into a baryon asymmetry through (hot) unsuppressed 
sphaleron transitions. If the PT is strongly first order this asymmetry is not 
destroyed by the Higgs phase (weak) sphaleron when finally the equilibrium phase 
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takes over . The rise of the generated baryon asymmetry depends on the spatial 
variation of the CP-violating phase in the wall. Moreover, in the MSSM the 
variation A/3 ~ d/3/dx of tan/3 is an important ingredient for the determination 
of the baryon asymmetry rj. Due to its smallness the observed baryon asymmetry 
may be reached only in a small, eventually tuned range [15] of parameters and 
therefore also a better knowledge of the wall velocity is needed, since there is 
a strong dependency on it. For instance in [16] it was found in semiclassical 
calculations that with a given Higgs field profile h 2 (x) 

[°° h'(x)d/}(x) 

"~L dx ^?— (L1) 

Thus the asymmetry increases with decreasing wall velocity. Of course v w — > 
is not possible, in this case there would be no out-of-equilibrium region. Quite 
naively one might expect that more particles would imply more interactions and 
in consequence a higher viscosity of the plasma. For this reason it usually was 
assumed (e.g. [17]) that the velocity of bubble walls of supersymmetric phase 
transitions is smaller than the velocity in the SM. But that has to be shown in a 
detailed calculation. 

To go beyond speculations we want to reconsider this question and the cal- 
culation of the wall velocity taking into account also supersymmetric particles. 
The proceeding is similar to that of [18]. Nevertheless, in case of two (Higgs-) 
scalars and arbitrarily many interacting particle species (fermions and bosons), 
the calculation and the results differ considerably from [18]. Therefore, we will 
explicitely outline the main steps of the calculation (see also [13]) to demonstrate 
the differences. On the other hand, at some points which can be found in [18] we 
keep short in the representation. 



2 Bubble Wall Equation of Motion 

Energy conservation leads to the equations of motion of an electroweak bubble 
wall interacting with a hot plasma of particles [18-20]: 

™+n*)+E^/pfb/.(M = o, (2.D 

where fa = fa,i + 5 fa is the distribution function for a particle species in the heat 
bath. We have to sum over all particle species i. The distribution function is 
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divided up into equilibrium part fo,i and out-of-equilibrium part 5 fa. For two 
scalars (or more, analogously) we obtain 

dV T (h 1 ,h 2 ) ^dm 2 f d 3 p 

dV T (h 1: h 2 ) ^dm 2 f d 3 p rr . . . . 

+ + ^mrj j^kw/f^ = °- < 2 - 3 > 

The equilibrium part has been absorbed into the equilibrium temperature de- 
pendent effective potential V^/ii,/^)- I n equilibrium we would obtain the free 
equations for critical bubbles respectively domain walls for large radii. 

In the following we will restrict ourselves to late times leading to a stationarily 
moving domain wall where the friction stops the bubble wall acceleration. This is 
the long period of bubble expansion where baryogenesis takes place. The influence 
of different friction or viscosity terms in Standard type models was investigated 
in various papers [19-21]. We assume not too large velocities and check the self 
consistency of this assumption in the MSSM. 



3 Fluid Equations 

Now we want to derive the deviations 5fi from the equilibrium population den- 
sities originating from a moving wall. We will therefore discuss the Boltzmann 
equation in the fluid frame, the "fluid equations": 

dtfi = dtfi + x^-fi + p x -£-fi = -C[fl (3.1) 
dx d Px 

with the population density /j and energy E = \fpl + m 2 (x). C[fi\ represents 
the scattering integral and will be discussed in section 6. The classical (WKB) 
approximation is valid for 

p > -?- ("thick wall"). (3.2) 

For particles with E,p>gT this should be fulfilled. Thus, here infrared particles 
are supposed not to contribute to the friction [18-20]. This is a crude approxi- 
mation. A further understanding of the infrared particle contribution is therefore 
needed which goes beyond the aim of this paper 3 . In the MSSM the wall thick- 
ness L w is of order 15/T-40/T, as found in [7,13,22], and L w ^> 1/T is fulfilled. 

3 Ref. [28] revisits the calculation of the friction in the SM. The prediction is that hard 
thermal loop effects play an important role in the damping of the gauge fields in the hot phase 
and that such coherent gauge field contributions are very effective for generating a friction. 
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With [18] we denote those particles which couple very weakly to the Higgs as 
"light particles" . Particles coupling strongly to the Higgs are heavy in the Higgs 
phase and therefore called "heavy" . However, "superheavy" particles as the "left 
handed" stops do not appear in our calculation besides their remnants in the 
effective potential. We treat as "heavy" particles only top quarks, (right handed) 
stops, and W bosons. Besides for the physical stop mass, we neglect the U(l) 
and treat only the SU(2) with its coupling g. The Higgses are left out. 

We assume now that the interaction between wall and particle plasma is the 
origin of small perturbations from equilibrium. We will treat perturbations in the 
temperature 5T, velocity 8v and chemical potential 8/j, and linearize the resulting 
fluid equations. Then the full population density of a particle species i in the 
fluid frame is given by 

where we have generally space dependent perturbations 5i from equilibrium. In 
principle one must include perturbations from the global value for each particle 
species. A simplification is to treat all the "light" particle species as one common 
background fluid. This background fluid obtains common perturbations 5v bg in 
the velocity and 5T bg in the temperature. This leads us to 



Si = - 



E 

5^ + —(STi + 5T bg ) + p x (5vi + 5v bg ) 



(3.4) 



for the "heavy" particles. The spatial profiles of all these perturbations depend 
on the microscopic physics. We treat particles and antiparticles as one species 
neglecting CP violation which is a minor effect on the friction. For the calculation 
of the baryon asymmetry this is the important effect and would be involved in 
eq. (3.4) by a perturbation 5E± in the energy dividing particles and antiparticles 
[4,23]. 

We now expand d t f to linear order in the perturbations. The Boltzmann 
equation can then be written as 



dtim 2 



(-/o) (^[9xSfi + ^d x (5T + 5T bg ) + Px d x (5v + 5v bg )} + d t 5fi 
+^d t (5T + 5T bg ) + Px d t (5v + 5v bg )^j + TC(5fi, 5T, Sv) = (-/£ 



2E 

(3.5) 

The term on the right hand side of (3.5) drives the population density away from 
equilibrium. 
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The collision term depends on all perturbations. But since the perturbations 
are Lagrangian multipliers for particle number, energy, and momentum, we can 
determine the parameters by the appropriate integration choice 

/d?p f d^p f d?p 

p' iVr and l v 'Wf- (3 ' 6) 

The resulting three equations coming from the Boltzmann equation, the "fluid 
equations", are coupled through the collision term C[5fi,5T,5v}. Performing the 
integrals leads to the general pattern 

/ wk^ c[f] = s ^ +STr ^ 

/ wk~^ c[f] = 5vTr - (3 - 7) 

where the rates T are of the form V ~ a 2 ln(l/a)T, and a = g 2 /(An) is the gauge 
coupling. 

The expressions (3.7) have to be evaluated graph by graph through the out of 
equilibrium interactions of each particle species. In sec. 6 we will calculate the 
leading contributions. 

For a stationary wall we can use d t Si — > f w <^, and d z 5i — » 5^, where the prime 
denotes the derivative with respect to z = x — v w t (Galilei transformation for 
small v w ). 

In the following the equations are again similar to those in [18] but there are 
important additional terms. We can solve the fluid equations to eliminate the 
background perturbations 5v' bg and ST^ : 

-v,(A + B) + C + D 
A + B-3vJC + D) 
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where we have used the notation 



A = N f (5fi f T, 2f + 5T f r T2f ), 

fermions / 
B = N b (5fx b T^ 2b + 5T b T T2b ), 

bosons b 
C = £ N f Sv,TT vf , 

fermions / 

D = N b$VbTT vb . (3.10) 

bosons b 

In contrast to the computations in [18] these terms may include arbitrarily many 
fermions and bosons, respectively. 

We can see that for walls moving with the velocity of sound v w — v s — 1/ ~ 
0.577, the approximation obviously breaks down. This demonstrates the limit of 
the expansion in linear perturbations. 

For each heavy particle species in the plasma we have three fluid equations 
resulting from the combination of eqs. (3.5), (3.6) and (3.7). The final form of 
the fluid equations can be written in a matrix notation: 

A5' + T5 = F, (3.11) 

where 

1 

r = r + — m. (3.12) 

c 4 

The matrices A, T, T , and M are given below. The number c 4 is the heat capacity 
of the plasma c 4 = 78c 4 _ + 37c 4+ including light quarks, leptons, and sleptons in 
the plasma. The number changes when we include further light particles. The 
perturbations are combined in a vector 5, the driving terms are combined in 
the vector F . The driving term containing (m 2 )' can be split up into different 
contributions 

The vector 5 and the matrices for k particle species (index x denotes + or — , for 
fermions and bosons, respectively, for the ith particle): 

5 =[5^ ST, T5 Vl ... SpL k 5T k T5v k ] T , (3.14) 
F = ^[c lx (mjy c 2x (mlY ... c Xx {m\y c 2x {m\)> of . (3.15) 
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The matrices A and T are of block diagonal form: 
fAi 



A 



r 



At 



ri 











C3i 

3 


where 


A, = 


v w c 3i 




Cii 

3 






CM 


C4i 








. 3 


3 


3 










" 


where 





















r 



(3.16) 



(3.17) 



Moreover, 



is a square matrix of the form 

"Mil ••• M u 



M M 



M 



fc,fe 



(3.18) 



Mjj = AT, 



(3.19) 



c 3ir M 2,j C3jr^ 2 ,j C3jr„j 
c 4i^fj,2,j C4iTx2,j C4iT v j 
_C4i^fj,2,j C4i^T2,j C4i^v,j_ 

where i in c^, denotes fermionic or bosonic contributions c%±, C4± of the ith and 
jih particle species, respectively. They are defined for fermions(+) and bosons(-) 
through 

E n ~ 2 d 3 p 



C n ± 



-/o(±): 



T"+ iJUV ^(2tt) 3 ' 

4 Higgs Equation of Motion and Friction 



(3.20) 



With the definition (3.20) and taking into account (righthanded) stop-, top- and 
W particles the equations of motion can be approximated in the fluid picture as 








in , 9V T/h , . N w dm 2 w 



_ 4 . + _£ (fcl ,fc l) + _ fr 



ci_5/i W + c 2 _<57V + c 2 -T5v 



^ dm? 
2 d/i 2 

_/V~ diTv? 

+ -ifT— A( Cl _^ + C2 _5T f + C2 _rH-). 



2 d/i 2 

This can formally be rewritten as 



(4.1) 



dVr T 

-K + ^(h l ,h 2 ,T) + -h 1 G l 5 l 

dVr T 

-ti2 + 7 n Z (hi,h 2 ,T) + -h 2 G2S2 

(7/7,2 Z 



o, 



(4.2) 
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with 



'N wCl .g 2 /2 
N w c 2 -g 2 /2 
N w c 2 .g 2 /2 




~N wCl -g 2 /2 
N w c 2 _g 2 /2 
N w c 2 „g 2 /2 
Nici-y 2 t 




5fi w 




5nw 






5T W 




5T W 






T5v w 




T5v w 



















, G 2 = 


N' t c 2 -y 2 t 




m 


, and 5 2 = 














T5v~ t 












N t c 1+ y 2 




5/it 












N t c 2+ y 2 t 




ST t 












N t c 2+ y 2 t 




_T8v t _ 








(4.3) 

The vectors 5±, 5 2 in eq. (4.3) contain the perturbation functions which can be 
found as solution of eq. (3.11). Here we used the splitting of the perturbation 
vector we discussed at eq. (3.13). Thus one is lead to a system of linear, coupled 
differential equations. It can be solved numerically. The solutions give the profiles 
of the perturbations in the wall. Such profiles are discussed in [21] for single scalar 
models. We now use again the thick wall limit. First we now approximate 5' = 0. 
(This will be improved later on). Then we can solve eq. (3.11) for the perturbation 
vectors 5i t2 : 

Si,2 = r~ 1 F li2 . (4.4) 
We obtain F 12 for stops, top and W from eq. (3.11) in the simplified form 



F 2 = ^ 
2T 



and 



2T 



ci-(m^)' 
c 2 -(m 2 v )' 





\ci-g 2 /2] 
c2-g 2 /2 




ci-(mf)' 
c 2 _(m|)' 



= ^h 2 h 2 


ci-y 2 
C2~y 2 




ci+(m?)' 




ci+y 2 


c 2+ (m 2 )' 





c2+y 2 




ci-(m^)' 
c 2 -(m 2 v )' 



Vw !, 1.1 

= — hiki 
2T 1 


~ci-g 2 /2 
c 2 -g 2 /2 




^h 2 h' 2 F 2 , 



(4.5) 



= -Ufa fo' p 
2T 1 



(4.6) 



Finally, we arrive at the bubble wall equations of motion with friction, 

h'l - V± 



h 2 



h! 2 ' - V± 



h 2 

T)2V w — h 2 . 



(4.7) 
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The dimensionless constants rji and r} 2 are denned through the relation 

Vi = ^G 1 T- 1 F 1 , (4.8) 

m = jG 2 T- 1 F 2 , (4.9) 

where we used the definitions for G lj2 in eq. (4.3) and F lj2 from eqs. (4.5) and 
(4.6) above. The factor T rescales ~ T _1 . The friction terms % and rj 2 
are still slightly tan/3-dependent. The constants give the viscosity of the medium 
which is perturbed by the moving wall surface. But since 8' in general is not 
negligible we need the full solution to eq. (3.11) and 77^2 cannot be defined as in 
eqs. (4.8) and (4.9) and the r.h.s. of (4.7) maintain an implicit t> w -dependence 
beside the factor. We then have to solve eq. (4.2) directly instead of eq. (4.7). 
In the next section we proceed with (4.7)in order to obtain an analytical formula 
for the wall velocity which is useful in any case. 

5 Wall Velocity in the MSSM 

In order to solve eqs. (4.7) we derive a virial theorem, based on the necessity that 
for a stationary wall the pressure to the wall surface is balanced by the friction: 
the pressure from inside the bubble which is responsible for the expansion and the 
pressure resulting from the viscosity of the plasma must be equal. The pressure 
on a free bubble wall can be obtained from l.h.s. of the equation of motion (4.7) 
by 

iK - gj^J Kdx = VrihiO)) - Vrihix = 00)) = AV T , (5.1) 
which is compensated due to the friction term, the r.h.s. of (4.7): 

/ (5 " h[dx = I ^ Vw Y {K)2 = AVt > (5 ' 2) 

where AVt is the difference in the effective potential values at the transition 
temperature T n , which is basically the nucleation temperature. 

Since we have two equations of the same type for each of the Higgs scalars and 
both of them develop friction terms, we have to add the pressure on the bubble 
surface. They are different due to the different particle species and couplings. 
Thus we find different pressures which would lead to different wall velocities, 
if the equations were completely decoupled. But due to the effective potential 
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which couples the equations of motion we have back-reaction. This leads to a 
change in A/3 = max(9 '(3/dz). It might be interesting to investigate this question 
in more detail, since a larger A/3 is highly welcome to obtain a larger baryon 
asymmetry. This point gets even more important with the knowledge of the re- 
sults of [13, 14] where we realized that in the MSSM transitional CP violation 
does not occur. Therefore we must exploit the explicit phases which may never- 
theless be strongly restricted by experimental bounds. The determination of A/3 
may be done numerically by solving eqs. (4.7) with extensions of the methods 
of [13,24,25]. 

For our estimate of the wall velocity we will use as approximation a constant 
tan/3 since the deviation is so strongly suppressed as found in [7,22,24]. Then we 
can add the expressions for the pressure for both Higgs fields, pi and p 2 : 

Pi+P2 = y{^ J h 2 1 (h' 1 ) 2 dx + r ]2 j h 2 2 (h' 2 ) 2 dx^ = 2AV T . (5.3) 



Next, we define an average direction and introduce a field h = \/h\ + h 2 . We 
find 

OAT/ T r 

T = sin 4 /3(77 2 + cot 4 /^i) / h 2 {h'fdx. (5.4) 
v w J 

We realize that the term corresponding to 771, in comparison to the second, is 
strongly suppressed. Already moderately small values of tan/3=2,(3,6) cause 
roughly a suppression of order 10,(100,1000). The wall velocity is therefore pre- 
dominantly determined by the second Higgs. Moreover, this behaviour also is 
supported by the strong Yukawa couplings of the stop and the top which couple 
asymmetric in favour of h 2 . A 'large' tan/3 > 2 leads to a friction term, which 
is solely determined by i] 2 . The limit tan/3 — > leads to a large but finite ve- 
locity determined by the (tan/3 dependent) rji. Nevertheless, due to a different 
particle content 771 does not recover SM friction in this limit. However, the lower 
experimental limit is tan/3 > 2. 

In [22, 24] it was realized that in the MSSM the approximation to the bubble 
wall profile by a kink is a rather good choice. Hence, we set for the common 
Higgs field h(x) 

/^) = ^(l +t anh^Y (5.5) 

where L w is the wall thickness and h crit is the nontrivial vev at the nucleation 
temperature T n . The integral in eq. (5.4) is evaluated as 

*W^ = y|J, (5-6) 
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leading to the desired equation for the wall velocity 

20L W AV T (T n )T n 
Vw ^sm^ + ^cot^)" 

The missing numbers for L w , h cr u, T n , and AV(T n ) can be independently deter- 
mined with methods described in [13,22,24,26]. 

6 Viscosity and Wall Velocity 

The next step is to determine the specific number for the friction terms. There- 
fore we have to calculate scattering and annihilation rates resulting from the 
corresponding Feynman graphs. First we consider strong interactions. These 
processes contribute to order g 4 . We will drop all terms of order m/T and there- 
fore also s-channel processes which are of order m 2 /T 2 (cf. [18]). The scattering 
amplitudes are given by the following graphs. In principle, a stop can scatter off 
quarks, squarks and gluons: 




The annihilation graphs can be obtained by exchanging the s and t direction. In 
our approximation, there are only a few contributions to the matrix elements. 
For instance, we assume all the squarks besides the stop to be superheavy and 
decoupled. Hence, the first graph in (6.1) would only contribute for the stop-stop 
scattering. But, as the stops are at the same temperature and velocity, there is no 
change in the perturbations: J^<5 = Ofi+(E p — E p >)5T+(p z — p' z )5v = 0. Therefore, 
stop-squark scattering does not contribute to the rates in this approximation. 

The corresponding types of graphs have to be taken into account for the rates 
of the top-quark. They can scatter off gluons, quarks and squarks. Compared 
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to the SM, a new scattering graph appears for the top which in the MSSM may 
scatter off squarks. 

We calculate the rates which are defined in eq. (3.7) by the integration of the 
collision integral C[f], hence we need 

C[/J - J {2^2E p 2E p> 2E k> ^ * 
x(2ir) 4 5(p + k-p' - k')V[f] 
with V[f] = f p f k (l±f p ,)(l±f k ,)-f p ,f k ,(l±f p )(l±f k ). (6.3) 

The sign depends on the statistics of incoming an outgoing particles. For stop 
gluon interactions we have a minus sign in any case. The integrals are calculated 
in the relativistic limit. In this approximation we find the following matrix el- 
ements (after averaging over outgoing particles). For annihilation of stops into 
gluons/ghosts we have 



137 4 4 st + t 2 
—g s +36g s — 



\Ma\ 2 = —gt + 36gt—^, (6.4) 



and for stops scattering off gluons/ghosts, 



,2 128 4 oo i st + S 



\M s \ 2 = —gt + S2gt—^- (6.5) 

which can be understood easily by crossing the annihilation result and taking into 
account that stops scatter off ghosts and anti-ghosts. The large numbers result 
from large color and momentum factors each. We include thermal masses of the 
exchange particle m p in the denominator to regularize the divergence for small t. 
Different from [18] we we fully integrate the rates numerically. Some details of 
the integration of the rates are given in App. A. For the numerical evaluation of 
the resulting integrals we use g 2 /(47r) = a s = 0.12, aw = 1/30. The integrals for 
top-W-scattering can be found in [18]. Nevertheless, we have to reevaluate those 
expressions with different plasma masses for tops, W bosons, and stops since the 
super symmetric plasma differs from the SM plasma [27]. We use as plasma mass 
of the stop for our scenario [12] 



ffl; 



2 1 1 

2 „,2 , „2 i 1 u 2 1 / a2 



- -< + (-gi + -hf - -(Af + ^)/m' Q )T\ (6.6) 

As a required condition for baryogenesis we need very negative —mfj. As plasma 
mass for the gluon we take m 2 g 2/3g^T 2 . We use m w = 5/6g 2 T 2 for the W- 
mass, ra\ = (l/Qg 2 s + 3/32g 2 + l/8h 2 )T 2 for the top mass [27]. Annihilations 
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Figure 1: Wall velocity in dependence on the parameter tan/3(T = 0) for itlq = 
2TeV, A t = /2 = 0, and m A = 400GeV for ml = -60 2 , 0, 60 2 GeV 2 . Lower bunch 
of graphs for 8' — 0, upper for 8' ^ 0. 

into leptons and quarks include lepton masses mi = 3/32g 2 T 2 and bosons with 



After eliminating the perturbations with help of the rates as shown above in 
eqs. (3.11), (4.4) and finally eqs. (4.8) and (4.9), we can determine rji and r] 2 . 
Numerically we find for tan/3 > 0(2) and m v < O(0 GeV) 



This emphasizes the domination of r] 2 . In the Standard Model with only one 
viscosity constant one finds roughly i] SM ss 3. 

The wall velocities are to be determined from the value of tan/3 at the transition 
temperature T n and the corresponding values of L w , h crit and the difference in 
the potential heights AVr(T n ) from the one-loop resummed potential. In the 
whole scenario we assume the left handed stops to be decoupled and the right 
handed ones very light. We have no mixing A t — /j, — 0, mq = 2TeV, and 
m A = 400GeV. For a strong PT we prefer small right handed stop parameters 
mfj. The value mfj = — (60GeV) 2 , thus m~ t = 161GeV, is very near to the triple 
point [5,26]. There we find a strong phase transition with v/T = 0.96 at one-loop 



mass m 2 = 
l/41nm 2 . 



l/Qg 2 T 2 leading to an effective plasma mass In <m 2 >= 3/4 In m 2 



+ 



m w O(0.1), 

T] 2 < 0(100). 



(6.7) 
(6.8) 
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level for tan/3 = 2.0. At two-loop level it is even larger permitting larger tan/3 for 
the same strength [5,6]. In this parameter range, eventually including some stop 
mixing, one can reach Higgs masses up to ran ~ 110 GeV which agree with the 
present experimental bound. 

Performing all the calculations the wall velocity turns out to be of the order 
of v w ~ 0.05 — 0.1 which is lower than in the Standard Model. Due to the large 
numerical factor of the matrix elements, the overall rates are larger than expected 
but probably sufficiently small for baryogenesis. In Fig. 1 the dependence on 
the parameter tan/3(T = 0) is shown in the physically interesting range 2 < 
tan/3 < 6 for two values of m v = — 60 2 , 0, 60 2 GeV 2 . The lower bunch of curves is 
determined for our first approximation 5' = the upper from the full numerical 
solution to eq. (3.11) with 5' ^ 0. We find 4 corrections of roughly a factor of three 
for small my (light stop). Very heavy stops should decouple more and more which 
would lead to to increasing wall velocities again. This behaviour can reproduced 
with the full solution. In Fig. 2 this principle behaviour is demonstrated over a 
wide range of m\j and for tan/3 = 3. Nevertheless, for large (positive) mjj 
the used approximations become worse and corrections of order 0{m 2 /T 2 ) are 
important. Then our approximations break down. However, this range is in any 
case physically disfavoured for a strong PT. 

7 Conclusion and Discussion 

The velocity of walls of expanding bubbles during the phase transition is an im- 
portant ingredient for the calculation of the baryon asymmetry at the electroweak 
scale. In this paper we developed the utilities for the calculation of the wall ve- 
locity in extensions of the Standard Model and applied them to the MSSM. We 
generalized the method of [18] to include more than one Higgs field and arbitrarily 
many particles and applied it to the MSSM with a light right handed stop. 

We estimated the effect of the stops on the wall velocity in the plasma of 
the simplest supersymmetric extension of the SM with all particles heavy or in 
equilibrium besides the light stop, the W bosons and the top quark. With our 
calculation we give estimates on the friction or viscosity coefficients of this plasma. 
There are two parts, one for the two CP even Higgs fields each. In the region 
of tan/3 which is interesting for baryogenesis the velocity depends moderately on 
tan/3 and we find values of v w « (5 — 10) • 10 -2 . We can answer the question of 

4 Since eq. (3.11) depends through A on v w we solved it iteratively. 
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Figure 2: Wall velocity for increasing stop mass parameter my. At very large my 
(physically disfavoured) the stops decouple from the plasma leading to a large 
velocity while the first approximation 5' = contrarily leads to decreasing vw 
The diagram is calculated for tan/3 = 3, A t — /j, — 0, and m A = 400GeV. 

the title with Yes: stops do slow down the bubble wall. 

We have not included effects of further out-of-equilibrium SUSY particles like 
charginos and neutralinos and gave only the leading order results. Gluinos may 
have considerable influence due to their multiple interactions. We estimated the 
possible effect of a light Higgsino which produce logarithmically dominating equi- 
libration rates. There may be a strong effect, reducing the bubble wall velocity. 

Also a more precise consideration of the mass of the outer legs in the graphs 
presumably give an effect. We calculated the matrix elements in the massless 
limit, but especially for the stop this might be questionable. Moreover we ne- 
glected possible effects on the transition parameters due to a space dependent 
background temperature. A change in the nucleation temperature T n and poten- 
tially large reheating might change the results. These effects may also affect a 
deviation of the bubble wall profile from the kink Ansatz, as discussed in [18]. 

The contribution of the stop in our WKB calculation appears to be at least 
of the same order of magnitude as the soft gauge boson contribution estimated 
in [28]. A further calculation using combined techniques including Bodeker's 
effective theory [29] were nevertheless useful. 

Low wall velocities agree with the requirements of baryogenesis and enlarge 



15 



the baryon asymmetry. Our results agree with the possibility that electroweak 
baryogenesis is a realistic scenario for baryon asymmetry of our Universe. 
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Appendix A Strong interactions 

The matrix elements for stop-gluon interactions after averaging over outgoing 
particles are 

I-Ma| 2 = ^ s 4 + 36^^!, and \M S \ 2 = ^fgi + ?>2gf-^ (A.l) 

for annihilation and scattering, respectively. These matrix elements have to be 
integrated for the rates. We use the formalism and the definitions of [30], App. 
A, for both, annihilation (s-channel) and scattering (t-channel) rates. 

Annihilation We rewrite s, t, u with the following identities 

s = u 2 — q 2 , u = —s — t, k = uj — p, k' — cu—p', (A. 2) 

t = ^ { [(p - k){p' - k') - q 2 ] + cos0 v / (4^-s)(4p'A;'-s)} . (A.3) 

Then the rates for annihilation are given by 

r M i = ^ j ' dR2\M A \ 2 {p',P,q,uA), (A.4) 

Tti = ^1 / ' dR{E p + E k )\M A \ 2 {jt,P,q,w,<l>), (A.5) 
I> = ^ [ JdR2E p \M A \ 2 (p / ,p,q,u; :( f ) ) : (A.6) 
T T2 = -L J dRE p (E p + E k )\M A \ 2 {p',p,q,u,<j)), (A.7) 
r„ = J dRp z (p z + k z )\M A \ 2 (p',p,q,u;,<l>), (A.8) 

where 

/i p2tt poo pu> f^i^ f'^T 

dR= Wj d<P j dU j dvJ^dPj^dp'fpMi + fp'^ + h')- 

(A.9) 
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Scattering Here we rewrite s, t, u with 



t = to 2 — q 2 u = —t — s, p' = p + tu, k' = k — cu, (A. 10) 

-t 
2q* 



s = ^ {[(p + p')(k + k') + q 2 ] - cos <f>y/(W + t)(4kk' +7)} .(A.ll) 



Then the rates for scattering are given by 

T T2 = -L J dRE p {E p -E p ,)\M s \ 2 {p',P,q,uA), (A.12) 

r« = 7^ J dRp z (p z -p' z )\Ms\\p\p,q^,<P) } (A.13) 

where 

dR= ¥^J d<p j, dq j d ^J w _ q d PJ w+q dk fpfk^+fp')^+fk')-(^) 



2 2 



The (ii? integration must be done numerically and leads to friction coefficients 
771 ~ 0(0.1) and r] 2 < (9(100), slightly depending on the physical parameters of 
the system. 
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